Search for: All records

Abstract Cosmic rays (CRs) interact with turbulent magnetic fields in the interstellar medium (ISM), generating nonthermal emission. After many decades of studies, the theoretical understanding of their diffusion in the ISM continues to pose a challenge. This study numerically explores a recent prediction termed “mirror diffusion” and its synergy with the traditional diffusion mechanism based on gyroresonant scattering. Our study combines 3D MHD simulations of star-forming regions with test particle simulations to analyze CR diffusion. We demonstrate the significance of mirror diffusion in CR diffusion parallel to the magnetic field when the mirroring condition is satisfied. Our results support the theoretical expectation that the resulting particle propagation arising from mirror diffusion in combination with much faster diffusion induced by gyroresonant scattering resembles a Levy-flight-like propagation. Our study highlights the necessity to reevaluate the diffusion coefficients traditionally adopted in the ISM based on gyroresonant scattering alone. For instance, our simulations imply a diffusion coefficient ∼10²⁷cm²s^–1for particles with a few hundred TeV within regions spanning a few parsecs around the source. This estimate is in agreement with gamma-ray observations, which show the relevance of our results for the understanding of gamma-ray emission in star-forming regions. 

Abstract The dynamics of star-forming gas can be affected by many physical processes, such as turbulence, gravity, supernova explosions, and magnetic fields. In this paper, we investigate several nearby star-forming regions (Orion, Upper Sco, Taurus, and Perseus) for kinematic imprints of these influences on the newly formed stars. Using Gaia DR3 astrometry and APOGEE DR17 radial velocities, we compute first-order velocity structure functions (VSFs) of young stars in galactic Cartesian coordinates in both 6D (3D positions and 3D velocities) and 4D (3D positions and each 1D velocity) to identify signatures of turbulence and anisotropic motion. We also construct 3D and 1D radial velocity profiles to identify coherent expansion trends, and compare stellar proper motions to plane-of-sky magnetic field orientations in Taurus and Perseus. We find that the VSFs are mildly anisotropic, with slightly different amplitudes, slopes, or features in different directions in several groups, but in general, they are all consistent with Larson’s Relation at intermediate length scales, especially in less compact groups. In several cases, the VSFs exhibit features suggestive of local energy injection from supernovae. Radial velocity profiles reveal clear anisotropic expansion in multiple groups, with the most extreme cases corresponding to those with the most anisotropic VSFs. In Perseus, we find that the motions of young stars are preferentially perpendicular to the local magnetic field. We find multiple, overlapping causes in each group for the observed kinematics. Our findings support that young stars remember more than just the turbulent state of their natal clouds. 

Abstract Magnetic fields and turbulence are fundamental to the evolutions of galaxies, yet their precise measurement and analysis present significant challenges. The recently developed Velocity Gradient Technique (VGT), which capitalizes on the anisotropy inherent in magnetohydrodynamic (MHD) turbulence, represents a new method for mapping magnetic fields in galaxies using spectroscopic observations. Most validations of VGT thus far have relied upon idealized MHD turbulence simulations, however, which lack the more complex dynamics found in galaxies and galaxy mergers. In this study, we scrutinize VGT using an AREPO-based cosmological galaxy merger simulation, testing its effectiveness across pre-merger, merging, and post-merger stages. We examine the underlying assumptions of VGT and probe the statistics of gas density, velocity, and magnetic fields over time. We find that the velocity fluctuations are indeed anisotropic at each stage, being larger in the direction perpendicular to the local magnetic field, as required by VGT. We find additionally that galaxy mergers substantially intensify the velocity and density fluctuations and amplify the magnetic fields at all scales. The observed scaling of the velocity fluctuations shows a steeper trend thanr^1/2between 0.6 and 3 kpc and a shallower trend at larger scales. The scaling of the magnetic field and density fluctuations at scales ≲1.0 kpc also predominantly aligns withr^1/2. Finally, we compare results from VGT to those derived from polarization-like mock magnetic field measurements, finding consistent and statistically significant global agreement in all cases. 

ABSTRACT We investigate the driving of MHD turbulence by gravitational contraction using simulations of an initially spherical, isothermal, magnetically supercritical molecular cloud core with transonic and trans-Alfvénic turbulence. We perform a Helmholtz decomposition of the velocity field, and investigate the evolution of its solenoidal and compressible parts, as well as of the velocity component along the gravitational acceleration vector, a proxy for the infall component of the velocity field. We find that (1) In spite of being supercritical, the core first contracts to a sheet perpendicular to the mean magnetic field, and the sheet itself collapses. (2) The solenoidal component of the turbulence remains at roughly its initial level throughout the simulation, while the compressible component increases continuously, implying that turbulence does not dissipate towards the centre of the core. (3) The distribution of simulation cells in the B–ρ plane occupies a wide triangular region at low densities, bounded below by the expected trend for fast MHD waves (B ∝ ρ, applicable for high-local Alfvénic Mach number MA) and above by the trend expected for slow waves (B ∼ constant, applicable for low local MA). At high densities, the distribution follows a single trend $$B \propto \rho ^{\gamma _{\rm eff}}$$, with 1/2 < γeff < 2/3, as expected for gravitational compression. (4) The mass-to-magnetic flux ratio λ increases with radius r due to the different scalings of the mass and magnetic flux with r. At a fixed radius, λ increases with time due to the accretion of material along field lines. (5) The solenoidal energy fraction is much smaller than the total turbulent component, indicating that the collapse drives the turbulence mainly compressibly, even in directions orthogonal to that of the collapse. 

ABSTRACT The coupling state between ions and neutrals in the interstellar medium plays a key role in the dynamics of magnetohydrodynamic (MHD) turbulence, but is challenging to study numerically. In this work, we investigate the damping of MHD turbulence in a partially ionized medium using 3D two-fluid (ions + neutrals) simulations generated with the athenak code. Specifically, we examine the velocity, density, and magnetic field statistics of the two-fluid MHD turbulence in different regimes of neutral-ion coupling. Our results demonstrate that when ions and neutrals are strongly coupled, the velocity statistics resemble those of single-fluid MHD turbulence. Both the velocity structures and kinetic energy spectra of ions and neutrals are similar, while their density structures can be significantly different. With an excess of small-scale sharp density fluctuations in ions, the density spectrum in ions is shallower than that of neutrals. When ions and neutrals are weakly coupled, the turbulence in ions is more severely damped due to the ion-neutral collisional friction than that in neutrals, resulting in a steep kinetic energy spectrum and density spectrum in ions compared to the Kolmogorov spectrum. We also find that the magnetic energy spectrum basically follows the shape of the kinetic energy spectrum of ions, irrespective of the coupling regime. In addition, we find large density fluctuations in ions and neutrals and thus spatially inhomogeneous ionization fractions. As a result, the neutral-ion decoupling and damping of MHD turbulence take place over a range of length-scales. 

Abstract The interstellar medium (ISM) is turbulent over vast scales and in various phases. In this paper, we study turbulence with different tracers in four nearby star-forming regions: Orion, Ophiuchus, Perseus, and Taurus. We combine the APOGEE-2 and Gaia surveys to obtain the full six-dimensional measurements of positions and velocities of young stars in these regions. The velocity structure functions (VSFs) of the stars show a universal scaling of turbulence. We also obtain Hαgas kinematics in these four regions from the Wisconsin H-Alpha Mapper. The VSFs of the Hαare more diverse compared to those of stars. In regions with recent supernova activities, they show characteristics of local energy injections and higher amplitudes compared to the VSFs of stars and of CO from the literature. Such difference in amplitude of the VSFs can be explained by the different energy and momentum transport from supernovae into different phases of the ISM, thus resulting in higher levels of turbulence in the warm ionized phase traced by Hα. In regions without recent supernova activities, the VSFs of young stars, Hα, and CO are generally consistent, indicating well-coupled turbulence between different phases. Within individual regions, the brighter parts of the Hαgas tend to have a higher level of turbulence than the low-emission parts. Our findings support a complex picture of the Milky Way ISM, where turbulence can be driven at different scales and inject energy unevenly into different phases.